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ABSTRACT 

The order perturbative QCD correction to the Gottfried sum rule is obtained. The 
result is based on numerical calculation of the order contribution to the coefficient 
function and on the new estimate of the three-loop anomalous dimension term. The 
correction found is negative and rather small. Therefore it does not affect the necessity 
to introduce flavour-asymmetry between u and d antiquarks for the description of NMC 
result for the Gottfried sum rule. 
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1 Introduction 



One of the still actively discussed problems of deep-inelastic scattering (DIS) is related to 
the consideration of the Gottfried sum rule P], namely 
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This sum rule is the /irst non-single (NS) moment of the difference of F2 structure functions 
(SFs) of charged lepton-nucleon DIS which in general has the following definition: 
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An extensive discussion of the current studies of this sum rule was given in the review of 
Ref. j2]. However, for the sake of completeness, we will remind the existing experimental 
situation,which is stimulating the continuation of the research of various subjects, related 
to the Gottfried sum rule. 

In fact if the sea is flavour symmetric, namely u=d, one should have 
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However, the most detailed analysis of muon-nucleon DIS data of NMC collaboration gives 
the following result 0: 



/^^^^(qs = 4 Qgv2) = 0.235 ± 0.026 
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It clearly indicates the violation of theoretical expression of Eq. and necessitates more 
detailed investigations of different effects, related to the Gottfried sum rule. In this Letter 
we reconsider the question of studying the contributions of a^, corrections to this sum rule 
previously raised in Ref. 



2 Available perturbative corrections 



The status of the 0(0;^) perturbative QCD corrections to Igsr "was summarised in Ref. 
[S]. Following this review, we will extend its presentation to the order level. 

It should be stressed that the renormalisation group equation for Iqsr contains the 
anomalous dimension term: 
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where As = as/ (47r) and 

/^^ = -2EA^:^' (6) 

The first two scheme-independent coefficients in Eq. (6) are well-known: 

/3o = (yC^-^n^) = 11 - 0.666667^^ (7) 
/?i = {^C\ - 2CFnf - ^CArif) = 102 - 12.6667n/ (8) 



3 ^ ' ' 3 
where Ci? = 4/3, Ca = 3 and is the number of active flavours. 

The corresponding anomalous dimension function has the canonical expansion 

7?=K = E7r'^:+' ■ (9) 

i>0 

However, like in the case of the first moments of SFs of uN DIS, the first coefficient of the 
NS anomalous dimension function of the first moment 7o^^ is identically equal to zero. 
The difference is starting to manifest itself from the two-loop level, where in order to get 
the corresponding result in case of anomalous dimension for Igsr it is necessary to make 
analytical continuation and to use the so-called (+) prescription (see e.g. Ref. 0). In the 
case of 7""^ this was done in Ref. ^ and Ref. jH] and result in the following analytical 
expression 

= -4(C^ - CfCa/2)[13 + 8C(3) - 271^] = +2.55755 (10) 

where the numerical value of C(3)=l. 2020569 was taken into account. The perturbative 
corrections to Igsr can be obtained from the solution of the renormalisation group equation 
ofEq. (5): 

Igsr{As) = ADiAs) X C{As) (11) 
where the anomalous dimension term is defined as 
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Since the first coefficient of 77jg'jj is identically zero ( namely 7q^^=0), there is no singular- 
ity in AD{As) and we can put in Eq. (fT^ the lower bound of integration 5 = 0. In this case 
we obtain the following expression for the expansion of AD{As) up to 0(a;^)-corrections: 

1 ^,11=1 1/1 ^'/^/"=l'^2 .71=1/3 . 71=1 \ 

A0M)-l^\\A,m.\{l^-2L^.^-^)Ai(Q^ . (13) 

The only unknown terms here is the third coefficient 72^^ of the anomalous dimension 
function 7y^^^(y4s), which in general is scheme-dependent. 

In the cases of n/ = 3 and n/ = 4 the numerical versions of Eq. (fT^ read 

AD{as)n,=, = l + 0.0355(^) + (- 0.0392 + |-^)(^)' (14) 
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AD{as)nf=4 = l + 0.0384(-^) + ( -0.0415 + ;^)(^) (15) 



where the scheme-dependent expression for 72 ^ is still unknown. Its value will be fixed 
in the next section using the results of calculations in the MS-scheme. 

Few words should be added here on the perturbative theory expansion of C{As). From 
the general grounds it should have the following form: 



l + Cr'As{Q') + Cr'Ai{Q 



(16) 



As was found in Ref. [9 its first coefficient is zero, namely C"~^ = 0. However, as will be 
shown in the next section the non-zero perturbative theory contribution is appearing at 
the two-loop level. 



3 Calculations and estimates of the a? contributions 



We will start from the calculations of perturbative contribution to the coefficient function 
C{As) at the a^-level. It can be obtained after applying (+) prescription to the results of 
Ref. [TU]. Indeed, the order correction to the coefficient function of Igsr is defined by 
taking the first moment from the sum 
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where the expressions for the functions C2^'*'*'~''(x, 1) and C2'^'^^\x^ 1) were calculated in 
Ref. jTU] and confirmed with the help of another technique in Ref. [11]. Integrating Eq. 
(fT7|) numerically with arbitrary Casimir operators Ca and Cp, we obtain the following 
n j-independent and scheme- independent result 
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Combining now Eq. (14) and Eq. (15) with Eq. (18) we find the following expressions for 
Igsr- 
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where ctg = as{Q'^) is the NLO expression for MS coupling constant. 

In order to get the feeling what might be the contribution of the terms proportional to 
72^^ we will avoid extrapolation procedure of the values of 73 used in Ref.[4j, calculated 
analytically for even n=2,4,..,14 in the works of Ref. P2]- Indeed, performing extrapolation 
from the even values of n for the NLO terms 7" of the corresponding anomalous dimension 
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function, we are obtaining the following estimate 7^^^ = 28.23, which is 10 times larger 
that the real value given in Eq. (10). Therefore, the used in Ref. jl] extrapolation procedure 
is considerably overestimating the value of the coefficient 7"^^. The similar situation can 
occur in the case of using extrapolation procedure for fixing the value of 72^^. Indeed, 
following the ideas of Ref. 4] we get from extrapolation of the known even values for 72 
the following estimates: 72^^ ?5i 361 for nj = 3 and 72^^ ~ 283 for Uf = 4, which to our 
point of view might be unrealistically large. 

Keeping in mind that only direct calculation of 72^^ can give the real numerical value 
of this term, we nevertheless are proposing the following way of fixing uncalculated contri- 
bution to 77^3'^ function. We noticed the following numerical pattern of the behaviour of 
anomalous dimension function for n >2: 7^/72 ~ 0.12 for ^/=4 (see Ref.|^ and Ref. |13j 
especially). We have checked that for nf=3 the similar relation is 7^/72 ~ 0.09 for n > 2. 
Hoping that these relations are also valid in case of n > 1 , we estimate the values for 
^^=1 = (1/0.12)77=^ « 21.3 in the case of n/ = 4 and 7^=^ = (1/0.09)7^=^ ^ 28.4 in the 
case of Uf = 3. Substituting them into Eq.(19) and Eq.(20) we get: 
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Taking now as(Q^) ~ 0.35 we arrive to the following numerical versions of Eq.(21) and 
Eq.(22): 
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Therefore in presented expression for the the order correction to the Gottfried sum rule 
is larger then the order a^-term. 

Theoretical errors to the presented third terms in Eqs.(21)-(34) are coming from the 
errors of 72^^ terms in Eqs. (19), (20), which is impossible to estimate without their 
direct theoretical calculations. In any case these terms are damped by huge numbers 
(64/5o) and it is unlikely that the direct calculations of 72^^ terms will change the results 
of Eqs. (23), (24) substantially. One can check this conclusion using the overestimated to 
our point of view results of application of the extrapolation procedure. Moreover, the 
main contributions to the a^-term in Eqs. (21)-(24) come from the calculated by us 
term of the coefficient function of the Gottfried sum rule. 



4 Comments on violation of the Gottfried sum rule 



In the previous section we found that order perturbative QCD corrections to the Got- 
tfried sum rule are really small and can not describe violation of the theoretical prediction 



4 



from its NMC experimental value. This, in turn, lead to the necessity of introduction of 
the effect of flavour asymmetry of antiquark distributions in the nucleon jS] , namely 
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This phenomeno logical result is important for fixing the corresponding d/u ratio in dif- 
ferent sets of parton distribution functions, which are relevant to the LHC physics (for a 
review see e.g. Ref. [15j ). On the other side the consideration of available E866 data 
for the Drell-Yan production in proton-proton and proton-deuteron scattering has con- 
firmed the effects of flavour asymmetry. Indeed the analysis of Ref. ^1] gave the following 
number 
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It was also noted in Ref. that it is unlikely to receive additional contribution to Eq. 
()26|) from the region above x = 0.35, since the sea is rather small in this region. However, 
the contribution to this whole integral from the unmeasured region x < 0.015 is missed. 
The attempt to fix it was made in Ref. ^H] using the extrapolation to small x region. As 
the result the authors of Ref. jTH] suggested the manifestation of substantial contribution 
of twist-4 1/Q^-effects in Eq. (1). Note, that final E866 result is 
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is closer to NMC result, than obtained in Ref.UHl extrapolated value, namely 
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0.09 ±0.02 



(2^ 



Therefore, in order to understand the status of their conclusion on the possibility of ex- 
istence of substantial contribution of the 1/Q^-corrections to the Gottfried sum rule it is 
necessary to be more careful in performing extrapolations to low x-region. It is highly 
desirable to estimate the effects of higher-twist contributions to the Gottfried sum rule 
using any concrete model. However, one should keep in mind that there are also some 
other explanations of the observed deviation from the canonical value 1/3 for the Gottfried 
sum rule (see Ref. and Ref. [2] for the review of other works on the subject). 



5 Conclusions 

In summary: we found non-zero O(a^) perturbative QCD contributions to the 
coefficient function of the Gottfried sum rule. We also estimated the effect due to non- 
zero value of the three-loop contribution to NS anomalous dimension function for n = 1 
moment, which as we think is rather small. More detailed result can be obtained after 
completing the analytical calculations of the three-loop corrections to the NS kernel of 
DGLAP equation. This work is now in progress (see e.g. Ref. JH])- In any case the 
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value of the ag-correction is dominated by the calculated in this Letter contribution to 
the coefficient function, which is negative, but also small. Therefore, the existing NMC 
observation of the flavour asymmetry between d and u antiquarks is surviving. We hope 
that the possible future new HERA data might be useful for more detailed measurement 
of the Gottfried sum rule and for further studies of the effect of flavour asymmetry in the 
d/u ratio. 
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